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ABSTRACT * 

We present a catalogue of 191 extragalactic sources detected by the Atacama Cosmology :■} 

Telescope (ACT) at 148 and/or 218 GHz in the 2008 Southern survey. Flux densities span 14 - y, 

1700 mJy, and we use source spectral indices derived using ACT-only data to divide our sources ® 

into two subpopulations: 167 radio galaxies powered by central active galactic nuclei (AGN) ^ 

and 24 dusty star-forming galaxies (DSFGs). We cross-identify 97 per cent of our sources (166 
of the AGN and 19 of the DSFGs) with those in currently available catalogues. When combined 
with flux densities from the Australia Telescope 20 GHz survey and follow-up observations 
with the Australia Telescope Compact Array, the synchrotron-dominated population is seen to 
exhibit a steepening of the slope of the spectral energy distribution from 20 to 148 GHz, with 
the trend continuing to 218 GHz. The ACT dust-dominated source population has a median 
spectral index, ai4g_2i8, of 3.7+o^g, and includes both local galaxies and sources with redshift 
around 6. Dusty sources with no counterpart in existing catalogues likely belong to a recently 
discovered subpopulation of DSFGs lensed by foreground galaxies or galaxy groups. 
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1 INTRODUCTION 

The technologies that enable observations of large numbers of mil¬ 
limetre and submillimetre sources were developed relatively re¬ 
cently. They open up for study a previously unexplored regime that 
has the power to reveal the evolution of underlying galaxy popula¬ 
tions over cosmic time, and in particular over epochs of intense star 
formation. Instruments such as the Submillimeter Common-User 
Bolometer Array (SCUBA, SCUBA-2; Holland et al. 1999, 2013) 
operating at 850 pm, the Balloon-borne Large Aperture Submil- 
limeter Telescope (BLAST; Pascale et al. 2008) operating at 250, 
350 and 500 pm, the Large Apex Bolometer Camera (Siringo et al. 

2009) operating at 870 pm and the AzTEC millimetre-wavelength 
camera (Wilson et al. 2008) operating at 1.1 and 2.1 mm have 
mapped up to 10 deg 2 of the sky, but greater coverage has been lim¬ 
ited by the large amount of integration time required to conduct blind 
surveys to significant cosmological depth. Increased sky coverage 
and sensitivity motivated the construction of space-based observato¬ 
ries such as Spitzer (Werner et al. 2004) and Herschel (Pilbratt et al. 

2010) , operating in the wavelength regime 3-500 pm and covering 
up to tens of square degrees. At longer wavelengths, the Wilkinson 
Microwave Anisotropy Probe ( WMAP\ Wright et al. 2009) covered 
the whole sky at frequencies up to 94 GHz (X = 3.2 mm), but 
observations of unresolved extragalactic sources were only com¬ 
plete above 2 Jy due to its large beam size (roughly 13 arcmin at 
94 GHz). The Planck space telescope contains a high-frequency in¬ 
strument observing at frequencies spanning 100-857 GHz (Lamarre 
et al. 2010). At 143 and 217 GHz (2.1 and 1.4 mm), the Planck 
beam sizes are approximately 7 and 5 arcmin, respectively, allow¬ 
ing for a source catalogue that is complete down to flux densities 
of about 1 Jy. Large-area (> 100 deg 2 ) ground-based radio surveys 
probe with higher resolution than WMAP or Planck , but only as 
high in frequency as 20 GHz (e.g. the Australia Telescope 20 GHz 
survey; Murphy et al. 2010). Thus, there is a niche for millimetre- 
wavelength large-scale mapping experiments with spatial resolution 
superior to that of the space-based observatories. 

One such experiment is the Atacama Cosmology Telescope 
(ACT; Swetz et al. 2011). The ACT collaboration released cosmic 
microwave background (CMB) temperature anisotropy maps made 
with arcminute resolution from its 2008 observing season at 148 
and 218 GHz (Diinner et al. 2013). These maps also contain galax¬ 
ies that are luminous at millimetre wavelengths. Measurements of 
the millimetre fluxes of a large sample of sources have great po¬ 
tential to discriminate among source population models and reveal 
new source populations. For example, recent millimetre surveys 
(e.g. Vieira et al. 2010; Marriage et al. 2011; Planck Collaboration 
VII 2013) have yielded source counts that led to updated source 
models (Tucci et al. 2011) and the discovery of a new class of high- 
flux dusty galaxies (this paper; Negrello et al. 2010; Vieira et al. 
2010). 

ACT-detected sources at flux densities greater than 20mJy are 
predominantly blazars that have synchrotron-dominated spectral 
energy distributions (SEDs). These galaxies are powered by active 
galactic nuclei (AGN), central black holes that accrete nearby matter 
in a process that produces time-variable jets out of the plane of the 
accretion disc; in the case of blazars, the jet direction lies close to our 
line of sight. Ejected ionized particles spiralling around magnetic 
field lines in the jets create the observed synchrotron emission. The 
study of blazar SEDs and source counts provides information about 
AGN physics (e.g. de Zotti et al. 2010; Planck Collaboration XV 
2011; Tucci et al. 2011). 


The second population identified at ACT wavelengths is com¬ 
prised of infrared-luminous, dusty star-forming galaxies (DSFGs), 
which exhibit modified blackbody emission at submillimetre to mil¬ 
limetre wavelengths, diminishing towards longer wavelengths. The 
observed SED is dominated by thermal emission from dust grains 
that have been heated by the prodigious optical and ultraviolet flux 
produced by newly formed stars (e.g. Draine 2003, and references 
therein). A fraction of these galaxies may have a significant non- 
thermal contribution from AGN at their cores, in addition to the 
thermal dust emission. 

Some of these DSFGs are local galaxies (z 1), cross-identified 
with Infrared Astronomical Satellite (IRAS) sources (Devereux & 
Young 1990). The IRAS wavebands (X = 12-100 qm) are relatively 
insensitive to dust with temperatures below 30 K, a significant and 
largely unexplored component of many nearby galaxies (Planck 
Collaboration XVI 2011). Recent results from Herschel (Amblard 
et al. 2010) and BLAST (Wiebe et al. 2009) have begun to extend our 
picture of the cold dust in galaxies, but millimetre-wavelength ex¬ 
periments can, through probing dust in nearby galaxies, contribute 
to establishing a well-calibrated dust SED for typical, low-redshift 
DSFGs. This in turn has cosmological ramifications, as current anal¬ 
yses of the first generations of stars and galaxies that fuel the cosmic 
infrared background (CIB) rely on understanding and extrapolating 
from dust SED templates. 

High-redshift DSFGs were observed by SCUBA in the first sys¬ 
tematic survey of these sources, which create a significant fraction 
of the CIB emission (Blain et al. 2002). Selected at 850 qm, these 
sources came to be known as submillimetre galaxies (SMGs). This 
new population of galaxies subsequently became a focus of ob¬ 
servations (e.g. Viero et al. 2009; Weifl et al. 2009; Austermann 
et al. 2010) and stacking analyses that resolved more of the CIB 
into emission from discrete, DSFGs (e.g. Dole et al. 2006; Devlin 
et al. 2009). Though most SMGs will be undetectable by current 
millimetre survey instruments, a new population of sources that 
are significantly brighter and rarer than the submillimetre-selected 
SMGs and that similarly exhibit dust-dominated spectral indices 
has been identified at millimetre wavelengths (Vieira et al. 2010). 
They do not have counterparts in the IRAS catalogue, indicating that 
they are not members of the standard, local, ultraluminous infrared 
galaxy (ULIRG) population. 

These sources have recently been shown to be a new, higher red- 
shift (z > 3) subpopulation of the progenitor galaxy background, 
brought to the fore in miIli metre-wavelength surveys because they 
are lensed by foreground galaxies or galaxy groups (Negrello et al. 
2007, 2010; Lima et al. 2010; Hezaveh et al. 2013; Vieira et al. 
2013). Such objects are extremely rare (at submillimetre wave¬ 
lengths, for example, see Rex et al. 2010; Lupu et al. 2012), but 
wide-area surveys will find more. The South Pole Telescope (SPT; 
Carlstrom et al. 2011) has reported significant numbers of these 
sources. This population provides an avenue for follow-up research 
to study the details not only of lensed SMGs, but also of the lens 
systems (e.g. Ikarashi et al. 2011 ; Scott et al. 2011 ; Lupu et al. 2012; 
WeiB et al. 2013). 

In this paper, we report for the first time the discovery of DSFGs 
in the ACT data. This is the first multifrequency analysis of the 
ACT sources and the second report on ACT extragalactic sources. 
Marriage et al. (2011, hereafter Mil) presented a catalogue of 
sources at 148 GHz only from the 455 deg 2 of the 2008 ACT South¬ 
ern survey with the best uniformity and coverage. Here we extend 
those results to include 218 GHz (and updated 148 GHz) flux den¬ 
sities for the same 2008 ACT Southern survey region from the new 
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(Diinner et al. 2013, hereafter D13) data release maps. Future work 
will present sources detected in the ACT equatorial strip. Mapmak¬ 
ing of the ACT 277 GHz data set is currently underway. 

The layout of the paper is as follows. Section 2 describes the 
2008 season ACT observations and the reduction of raw data into 
maps, as well as follow-up observations made with the Australia 
Telescope Compact Array (ATCA). Section 3 details our method 
of source extraction. The source catalogue, including its astromet¬ 
ric and flux density accuracy, and its estimated completeness and 
purity, is discussed in Section 4. Section 5 compares our catalogue 
with currently available data sets; Section 6 gives the source num¬ 
ber counts. Trends observed in the spectral indices of our source 
populations are analysed in Section 7. We conclude in Section 8 
with a summary of our results. 

2 OBSERVATIONS AND DATA 
2.1 ACT observations 

The ACT experiment (Swetz et al. 2011) is situated on the slopes 
of Cerro Toco in the Atacama Desert of Chile 1 at an elevation 
of 5190 m. ACT’s latitude gives access to both the northern and 
southern celestial hemispheres. Observations occurred simultane¬ 
ously in three frequency bands, at 148 GHz (2.0 mm), 218 GHz 
(1.4 mm) and 277 GHz (1.1 mm) with angular resolutions of roughly 
1.4,1.0 and 0.9 arcmin, respectively. Observations of Saturn were 
used to determine beam profiles and pointing (Hincks et al. 2010; 
Hasselfield et al. 2013). From 2007 to 2010, ACT targeted two sur¬ 
vey regions: the southern strip centred around <5 = —52?5 and the 
equatorial strip centred around <5 = 0°. Further information about 
the ACT observations can be found in D13. 


2.2 ACT data 

The reduction of raw ACT data into maps is detailed in D13; we will 
briefly review that process here. Each ACT detector array (one per 
frequency band) is composed of 1024 detectors, with each detector 
time stream first analysed and then kept or rejected based on multiple 
criteria, such as telescope operation, weather conditions, cosmic ray 
hits or other interference. Approximately 800 (700) h of data from 
the 2008 southern strip for 148 GHz (218 GHz) remain after these 

Maximum likelihood maps of pixels 30 arcsec on a side are 
produced from the time stream data. From an initial estimate of the 
maps, source profiles are fitted for. We then subtract the time stream 
models for the point sources from the data and re-map. This prevents 
point source power from being abased into the map, and improves 
the final flux density estimates for S/N > 5 sources (quantified be¬ 
low). Similarly, an estimate of the CMB signal is also removed. 
The mapping equation (Tegmark 1997) is then solved iteratively 
using the preconditioned conjugate gradient method. Finally, the 
signals that had been subtracted are added back to the map. We 
conducted simulations in which a signal corresponding to sources 
of known flux were injected into raw ACT detector time streams 
before mapping. The new mapping procedure with ‘source subtrac¬ 
tion’ described briefly here and more thoroughly in D13 (Section 
11.4) corrects a 3-5 percent downward bias in the recovered flux 
densities. 


1 22?9586 south latitude, 67?7875 west longitude. 


The measured source flux densities have converged on a sin¬ 
gle value by the 25th iteration of the maps used in our study. For 
148 GHz, the D13 release map is from iteration 1000 of the map- 
maker, and we use this map both for convenience and in order that 
this study be exactly reproducible using public data. For our study 
of source flux densities at 218 GHz, we use iteration 200, which is 
more than adequate. Between the map iterations used in this study, 
the fluxes of the sources change less than 1 per cent. One may ask 
whether the value to which the measured source flux density has 
converged is accurate. Based on end-to-end simulations in which 
mock sources are injected into the ACT time streams, we estimate 
this accuracy at 3 per cent. 

The 148 GHz map is calibrated in temperature through cross¬ 
comparison of spectra over the range 400 < l < 1000 with WMAP 
as in Hajian et al. (2011) with an uncertainty of 2 percent. This 
calibration is transferred to the 218GHz map (2.4 percent uncer¬ 
tainty) through cross-calibration with 148 GHz. Hasselfield et al. 
(2013, hereafter H13) used this calibration in measurements of the 
temperatures of Uranus and Saturn. The resulting planet temper¬ 
ature estimates were consistent with previous analyses (Griffin & 
Orton 1993; Goldin et al. 1997). To summarize the results of H13, 
both previous analyses adopted a flux standard based on a model of 
Mars (Wright 1976; Ulich 1981) with 5 percent systematic errors. 
The H13 estimate for the temperature of Uranus and Saturn was 
5 per cent below those of the previous analyses, a deficit which was 
also seen in the WMAP analysis of Weiland et al. (2011) at 94 GHz 
and is most likely due to systematic errors in the Mars standard. 
The agreement between the H13 analysis, calibrated on extended 
CMB emission, and the previous analyses, calibrated on Mars, lends 
confidence that the ACT WMAP-based calibration at large angular 
scales will also apply to the calibration of sources at small angular 

Photometry based on matched filtering (Section 3) relies on an 
accurate estimate of the instrument beam shape. The impact of beam 
error on photometry depends on the details of the filter (equation 1), 
which in turn depend on both the beam window function and the an¬ 
gular power spectrum of the background and noise. We have found 
that the error in the solid angle of the beam is a conservative proxy 
for the photometric error due to beam uncertainty, and we use this 
simpler quantity (solid angle uncertainty) in the beam-related pho¬ 
tometry error estimates described here. Due to uncertainties in the 
beam measurements from Saturn, we assign 1 percent photometric 
errors for both bands. The profiles of bright point sources in the sur¬ 
vey maps suggest that the effective beam for the survey is slightly 
broader than that measured from Saturn. The survey maps consist 
of overlapping observations taken over the entire season, and the 
observed broadening may be attributed to a night-by-night jitter 
in the telescope pointing with an rms of 5 ± 1 arcsec (H13). The 
uncertainty in jitter correction corresponds to 1 percent (148 GHz) 
and 2 percent (218 GHz) photometric uncertainties, which are cor¬ 
related between bands. 

Finally, for our photometric uncertainty budget, we account for 
the fact that AGN have a lower effective frequency band centre (and 
thus slightly broader beam) than Saturn, and DSFGs have a higher 
effective frequency band centre (and thus narrower beam; Swetz 
et al. 2011 , table 4). We choose effective 148 and 218 GHz frequency 
centres corresponding to half-way between a steep-spectrum AGN 
and a DSFG: 148.65 and 218.6 GHz. This choice introduces a pho¬ 
tometric bias of less than 1.5 percent at 148 GHz and 1.1 percent 
at 218 GHz, which is positive for DSFGs and negative for steep- 
spectrum AGN. We add in quadrature this photometric bias from the 
source spectrum to uncertainties due to mapping (3 per cent), WMAP 
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Figure 1. Sensitivity maps with source detections at 148 GHz (top) and 218 GHz (bottom), showing the most uniform 455 deg 2 patch of the southern strip 5- 

with the greatest depth of coverage. It lies between right ascension 00 h 12 m and 07 h 07 m , and declination —56° 11' and —49°00'. The deepest data correspond E. 

to an exposure time of 23.5 min arcmin -2 and a la sensitivity of 2.34 mjy at 148 GHz and 3.66 mjy at 218 GHz. White circles or squares mark the locations 
of ACT sources with a size proportional to the log of the associated source flux density. Circles denote sources designated as AGN, and squares denote sources ° 

with spectra indicative of DSFGs. Towards the edges of the map, the variation in local noise properties due to uneven coverage is more apparent. a- 


calibration (2 percent, 2.4 percent) and beam shape (1.4 percent, 

2.2 percent) to obtain an overall flux density calibration uncer¬ 
tainty of 4.1 percent at 148GHz and 4.6 percent at 218GHz. As 
shown in Section 5, the ACT flux densities agree with independent 
measurements to within this margin of error. 

For this study, we have used the most uniform 455 deg 2 of the 
2008 ACT southern strip at 148 and 218 GHz. This region, shown 
in Fig. 1, spans declination —56?2 < 5 < —49.0 and right ascen¬ 
sion 00 h 12 m < a < 07 h 07 m . At 148 GHz, we use the data publicly 
released with D13. 2 Typical white noise levels in this region of the 
map are 30 pK arcmin at 148 GHz and 50 pK arcmin for 218 GHz. 
As described in Section 3, when matched-filtered with the ACT 
beam, this white noise level results in a ltr point source flux density 
sensitivity in the best covered regions of <to = 2.34 mjy at 148 GHz 
and3.66mJyat218GHz. The sensitivity levels in Fig. 1 are propor¬ 
tional to the square root of the number of observations at that map 
location, N a b s , with one observation per 0.005 s. Then the sensitivity 
level in a given portion of the map is a = cr 0 y/ A^max/Nobs • 

2.3 ATCA observations 

Mil found that the sample of ACT 148-GHz -detected sources 
cross-identified with the Australia Telescope 20 GHz (AT20G) sur¬ 
vey (Murphy et al. 2010) is dominated by sources with peaked or 
falling SEDs using flux densities measured at 5, 20 and 148 GHz. 
The study also confirmed the findings of the AT20G study (Murphy 
et al. 2010; Massardi et al. 2011b), namely that this population of 
radio sources is characterized, on average, by spectral steepening 
between 20-30 GHz and 148 GHz. However, this sample from Mil 
was incomplete, biased in a way that favoured sources with negative 
spectral indices between 20 and 148 GHz due to the AT20G survey 
completeness level of 78 percent above 50mJy at 20GHz. ACT- 
selected sources with 148 GHz flux densities less than 50 mjy and 
flat or rising spectra may not have been detected by AT20G. There¬ 
fore, in order to complete the Mil 20-148GHz spectral study, a 
targeted set of measurements of flux densities at 20 GHz was made 


for the Ml 1 sources that were not identified with any source in the 
AT20G catalogue within a 30 arcsec search radius. 

Scheduling and weather constraints permitted us to observe 41 
of these 48 sources with the 6 x 22 m antenna array of ATCA 3 over 
6.5 h at 20 GHz on 2010 November 10, when the array was in its 
east-west 750A configuration, with 15 baselines ranging from 77 
to 3750 m. The primary beam of each array telescope at 20 GHz 
is 2.3 arcmin, with a resolution of 0.5 arcsec possible with the full 
array in this configuration. The synchrotron-dominated spectra of 
these sources, as revealed in M11, indicated that they were compact 
(AGN) and would be unresolved. For these observations, the two 
2-GHz-wide frequency bands of the new ATCA Compact Array 
Broadband Backend (Wilson et al. 2011) digital array correlators 
were set to be adjacent by centring them at 19 and 21 GHz. The 
reduced average flux densities over the whole bandwidth of the cor¬ 
relator correspond to the ‘20 GHz’ flux density. Good weather con¬ 
ditions prevailed throughout the observations, from 07:35:24.9 ut 
through 14:11:54.9 ut. A few data blocks were flagged and removed 
in order to minimize noise and any spurious effects. 

The primary flux calibrator used for these observations was 
PKS B1934—638, and the bandpass calibrator used was PKS 
B1921—293. The target sources were expected to have flux den¬ 
sities at or just below the AT20G survey limit of 40 mjy. Each 
target source was observed once for 1.5 min for an rms noise level 
of < 0.15 mjy beam -1 . Target observations were interleaved with 
observations of four secondary calibrators for pointing and phase 
corrections, chosen to lie close to the targets. 


2.4 ATCA data 

The ATCA follow-up data were reduced using a fully automated, 
custom, shell script pipeline based on tasks from the miriad aper¬ 
ture synthesis reduction package (Sault, Teuben & Wright 1995). 
An initial inspection of the data was performed to identify contam¬ 
ination or any problems in the data acquisition. Automatic proce¬ 
dures were used to identify and flag data for each frequency band 
affected by shadowing or known radio contamination, resulting in 
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2 http://lambda.gsfc.nasa.gov/product/act/ 


3 http://www.narrabri.atnf.csiro.au/ 
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less than 1 per cent of the band being flagged and cut. The pipeline 
then generated the calibration solutions for bandpass, flux density 
amplitude and phase based on the calibrators. After checking that re¬ 
sults were consistent between the two frequency bands, we merged 
the calibrated target visibilities before extracting the flux densities 
to improve sensitivity. 

Analyses of the AT20G survey (Murphy et al. 2010) and of other 
ATCA projects (Bonavera et al. 2011; Massardi et al. 201 la) have 
shown that the triple-correlation technique (Thompson, Moran & 
Swenson 1986) can be effectively used to determine source flux 
densities down to tenth of mjy scales for point-like sources unaf¬ 
fected by poor phase stability. This technique is effective for point 
sources even if they are not in the phase centre and in the case of 
low S/N. 

We manually inspected the 20 GHz deconvolved images, and 
wherever a source was clearly identifiable we measured its flux 
density and flux density error at its peak in the image and the rms 
of the image pixels in a region unaffected by the source emission. 
For sources with high S/N, deconvolution techniques were able to 
reconstruct images, despite the poor uv-coverage of our observa¬ 
tions. In the few cases where the source appeared extended with 
respect to the synthesized beam, we estimated the flux density by 
integrating over the source area. In some cases, the image showed 
a confused field with multiple peaks, primarily due to sidelobes; 
for these cases we extracted a flux density estimate using the triple¬ 
correlation technique. The flux density errors were obtained from 
the rms of the visibilities of the V Stokes parameter. This assumes 
that these objects have negligible circular polarization, which is an 
overestimate in the cases where the source is not in the phase cen¬ 
tre. The same statistic (V Stokes rms) was used to indicate the noise 
level reached for fields where no source was identified at 20 GHz. 
Of the 41 ACT sources observed in the follow-up campaign de¬ 
scribed here, all had a measured flux density at 20 GHz. The final 
results of this study are included as part of Table 4 with asterisks 
indicating sources that were part of the follow-up ATCA sample. 

3 ACT SOURCE EXTRACTION 

We use the same matched filtering method as described in Mil to 
produce an estimate of the amplitude of a source at the location of 
the source centre. In this section, we summarize the method with 
an emphasis on aspects for which this analysis differs from Ml 1. 

Initially, the data are weighted by the square root of the number 
of observations per pixel. This method results in a map with an 
approximately constant white noise level and a natural apodization 
as the number of observations falls off towards the edge of the survey 
region. Using a new procedure, we filter the entire D13 148 GHz 
release map and corresponding 218 GHz map, both of which extend 
beyond the deep 455 deg" area of this study, in order to eliminate 
any potential edge effects. We then use only the 455 deg" region for 
source identification. 

The following filter is applied to the Fourier transform of the 

*(*) =_ _ (1) 

/ B*(kf)F^ klt (kf)\ 7’ 0thcr (ft') I ~ 2 B(k') d k 

where k — (k x , k y ) is the angular wavenumber, and x and y re¬ 
fer to the right ascension and declination directions. B(k) is the 
Fourier transform of the ‘effective’ instrument beam. As described 
in Section 2, the ‘instantaneous’ beam is derived from planet ob¬ 
servations (H13). The ‘effective beam’ in the 2008 survey map is 
broadened by imperfect telescope repointing with mean deviation 



Figure 2. Matched filters (solid lines) are band pass filters that smooth 
the map on the scale of the beam (dashed lines) and remove large-scale 
structure associated with the CMB, other astrophysical signals (e.g. the 
Sunyaev-Zel’dovich effect) and residual contamination from atmospheric 
brightness fluctuations. The beam and filter functions are plotted with unit 
normalization at their peak (0 = 0). The matched filters were binned in 
radius to better show the relevant angular scales. In the analysis, we use 
two-dimensional filters to capture the anisotropic character of the noise. 


<j0 = 5 arcsec. This broadening is included in B(k) by multiplying 
the instantaneous beam transform (released with D13 for 148 GHz) 
by exp(— l 2 Og /2). T othC r is the Fourier transform of all components 
of the data besides point sources (i.e. atmospheric or detector noise, 
CMB, etc.). The function Fq, ^ (k) is a high-pass filter that removes 
undersampled large-scale modes below ko = 1000 and modes with 
| k x | < 100, which are contaminated in a fraction of ACT data 
by telescope scan-synchronous noise. While B(k) is well approxi¬ 
mated as azimuthally symmetric, we retain the full two-dimensional 
power spectrum | T„ lhcr (k') | 2 to downweight anisotropic noise in the 
maps. The azimuthally binned 148 and 218 GHz real space filters 
and associated beam profiles are shown in Fig. 2. 

The power spectrum | fotherC^O I 2 used in this analysis was con¬ 
structed in a different manner from Mil: we used the power spec¬ 
trum of the data itself instead of the average of difference (noise) 
maps and models for the CMB and other contaminating sky emis¬ 
sion. This estimate is robust since the total power from the extra- 
galactic source signal is low compared to the CMB, atmospheric 
noise and white noise. To avoid a noisy estimate of the power 
spectrum, we smooth the power spectrum | T athcr (k') | 2 with a Gaus¬ 
sian. The exact formulation of the smoothing does not significantly 
change the resulting filtered map. 

Applying the matched filter can cause ringing in the maps 
around the very brightest sources, which impacts source extrac¬ 
tion around other, low-S/N, sources. Therefore, we identify sources 
with S/N > 50 in an initial application of the filter, and mask them 
in the maps. Sources with S/N < 50 are then extracted by matched 
filtering these masked maps. Once filtered, groups of map pixels 
with S/N > 4.8 are identified as candidate sources, as illustrated in 
Fig. 3. 

Finite map pixel size affects the measured flux densities, since 
it is rare that a source falls exactly in the centre of a pixel, leading 
to a systematic negative bias and increased scatter. We increase the 
map pixel resolution by a factor of 16 in a region 0103 on a side 
around each source through Fourier space zero-padding (e.g. Press 
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Figure 3. Filtered 148 GHz (top) and 218 GHz (bottom) submaps. The data 
have been matched-filtered such that the grey-scale is in units of flux density 
(mjy) with white (black) corresponding to — lOmJy (30mjy). Insets show 
the flux density distribution across the data as a grey histogram that has a 
standard deviation of 2.57 mjy (3.78 mjy) for 148 GHz (218 GHz). Several 
sources, both synchrotron and dust dominated, marked as circle and square 
outlines, respectively, are apparent as black beam-sized flux excesses. The 
white extended object in the 148 GHz map, marked by a diamond, shows 
the Sunyaev-Zel’dovich effect decrement from Abell 3128 NE (ACT-CL 
J0330—5228). 


et al. 2007) to allow for a more precise determination of the source 
peak location, and therefore its flux density. Properly centring the 
detection has the effect of boosting the S/N of typical S/N = 4.8 
sources to >5. Lastly, to account for the convolution of sources 
with the map pixel, which acts as a low-pass filter (the square pixels 
become a sine function in Fourier space), we deconvolve the map 
pixel window function in the higher resolution map. We find that 
this method reduces systematic errors associated with pixelization 
to below 1 percent. Purity tests (Section 4.3) reveal a significant 
number of false detections below S/N of 5.0, due to local noise and 
striping artefacts. Therefore, we impose an S/N = 5.0 threshold 
on sources detected in each map assuming no prior knowledge of 
source location from the other frequency. 

Once an S/N >5.0 catalogue has been generated for each of 
the filtered 148 and 218 GHz maps, we look for sources that have 
been identified at one frequency but not the other. The flux density 
for the source at the second frequency is then measured from the 
filtered map for the second frequency. We use the same prescription 
as for the flux measurement method just described, using a submap 
at the second frequency centred on the location of the detection as 
determined from the map at the first frequency. 


4 THE 148 AND 218 GHZ CATALOGUE 

The ACT-detected source catalogue is given in Table 4. We find 169 
sources selected at 148 GHz with S/N > 5, spanning two decades in 
flux density, from 14 to 1700 mjy. The 218 GHz map independently 
yielded 133 sources with S/N > 5. The combination of these two 
independent source lists gives a total count of 191, with 110 galaxies 
detected with S/N > 5 at both frequencies. 

The catalogue provides the IAU name, celestial coordinates 
(J2000), S/N and flux density estimation of each 148 and 218 GHz 
ACT-detected source. Raw flux densities are estimated directly 
from the map as described in Section 3. Deboosted flux densi¬ 
ties as derived according to Section 4.2 are given with associated 


68 percent confidence intervals. The raw and deboosted spectral 
indices, ai48-2is, between 148 and 218 GHz for each source, are 
provided. The classification of sources as dust or synchrotron dom¬ 
inated is based on the spectral index criterion described 

in Section 4.2. If the source was cross-identified with an AT20G 
catalogue source (see Section 5), the AT20G source ID is given. If, 
instead, the 20 GHz source flux density was measured during the 
2010 November follow-up campaign (Section 2.3), the 20 GHz ID 
name has an asterisk next to it. Sources not cross-identified with 
one of the catalogues fisted in Section 5 are marked with a ‘d’ 
superscript. 

Correlating 148 GHz flux densities between this catalogue and 
the catalogue given in Mil, we find an average agreement at the 
2 percent level, with larger scatter for individual sources. This level 
of consistency is expected given changes in calibration, mapmaking, 
beam profile estimates and deboosting procedure made for this 
updated study. 

The reported flux densities are the average over approximately 
two months of observation for each source, many of which are 
AGN-driven radio galaxies and thus likely to have varied in that 
time. For example, of the three bright sources cross-identified with 
Planck , in the year between ACT and Planck observations one 
source varied by 20 percent in flux density whereas one did not 
vary at all within errors. However, our simultaneous multifrequency 
observations allow for a consistent internal spectral characterization 
between ACT bands. In a future study, multiple years of data will 
be used to quantify the effects of variability on individual source 
spectra. 

The following sections provide details and context for the cata¬ 
logue values. 

4.1 Astrometric accuracy 

Radio interferometers can achieve very precise positional accuracy 
for sources, so ACT-selected sources cross-identified with a robust 
radio catalogue give a good measure of the positional accuracy 
of the ACT source detections. The AT20G catalogue covers the 
southern sky, and through pointing checks against very long baseline 
interferometer measurements of International Celestial Reference 
Frame calibrators, the positional accuracy of AT20G is shown to be 
accurate to better than 1 arcsec (Murphy et al. 2010). We exclude 
nearby extended/resolved sources from this analysis, determined 
from cross-identification of our sources with currently available 
catalogues, as distant point-like sources will present a clearer picture 
of the overall accuracy of our pointing. 

We compared the positions of the 34 ACT sources with 
S/N148 > 16 to positions of associated sources in the AT20G cata¬ 
logue using a search radius of half the ACT 148 GHz beam (0.7 ar- 
cmin). Fig. 4 shows the offsets in location between the AT20G right 
ascension (RA) and declination (Dec.) and the ACT-derived posi¬ 
tions for the sources as they were detected in the 148 GHz (black 
points) and 218 GHz (blue points) maps. The results are also sum¬ 
marized in Table 1. The large overlaid error bars are centred on the 
mean offsets, and extend as far as the rms of the offsets. For sources 
with lower S/N, the ACT location rms with respect to the AT20G 
position becomes inflated by the effect of noise in the maps. There¬ 
fore, the catalogue uses coordinates derived from the 148 GHz map 
since there is a higher level of noise present in the 218 GHz map. 

4.2 Flux density accuracy 

Flux densities derived from the ACT maps have systematic uncer¬ 
tainties arising from five effects: the overall calibration uncertainty, 
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Figure 4. Astrometric accuracy of the ACT source detections. The 34 filled 
black circles denote the positional offsets of ACT sources with S/Nus >16 
from counterparts in the AT20G catalogue as calculated from the 148 GHz 
map. The turquoise triangles are for the same sources, but using positions 
derived from the 218 GHz map, which has more noise. The large error bars 
are centred on the mean offsets, and extend as far as the rms of the offsets. 
See Table 1. 

Table 1. Astrometric pointing accuracy. 



-10 -50 5 10 

Right Ascension (") 


148 GHz 218 GHz 


Mean RA offset (arcsec) 
Mean Dec. offset (arcsec) 
RMS in RA offset (arcsec) 
RMS in Dec. offset (arcsec) 


0.6 ± 0.4 0.5 ± 0.6 

-0.4 ± 0.3 0.0 ± 0.6 

2.1 3.5 


the mapmaker, errors in the assumed source profile, error in the 
assumed source spectrum and flux boosting of lower significance 
candidates. The first four potential sources of flux density error were 
discussed in Section 2. Here we will discuss the last source of flux 
density error. 

The differential counts of the sources selected in our sample fall 
steeply with increasing flux (Table 2, Fig. 6). With no prior infor¬ 
mation about the source flux, the most likely scenario is that the 
measured flux is the sum of a dimmer intrinsic flux and a positive 


noise fluctuation. We use the two-band Bayesian method developed 
in Crawford et al. (2010), and report the 16, 50 and 84 percentiles 
(68 percent enclosed, equivalent to lcr) of the posterior flux and 
spectral index distributions. For the source count priors in this cal¬ 
culation, we use the sum of the models of de Zotti et al. (2010) 
for radio sources (using the Tucci et al. 2011 model results in 
flux differences of <0.03<x) and Bethermin et al. (2011) for dust- 
dominated sources. Following Vieira et al. (2010), we take a flat 
prior on the spectral index between —3 and 5, consistent with the ex¬ 
pected range for our populations. The two-band likelihood includes 
negligible correlation between bands and is consistent with back¬ 
ground astrophysical emission rather than correlated atmospheric 
emission. 

The source populations in 148 and 218 GHz naturally split into 
sources having their emission dominated by synchrotron (centred 
on a = — 0.6) or thermal dust (centred on a = 3.7; see Fig. 7, bottom 
panel). We use the threshold spectral index a = 1.66 (Vieira et al. 
2010 ) to divide these populations in terms of their posterior spectral 
index populations, with P(a > 1.66) > 0.5 classified as dusty and 
P(a > 1.66) < 0.5 as synchrotron dominated. The classification is 
robust within a ± 0.5 of this threshold. While calibration uncertainty 
is included in all quoted fluxes, we conservatively ignore inter¬ 
band correlation of calibration errors when flux deboosting to avoid 
making assumptions about the correlation of the pointing jitter, 
chromatic effects and mapping/flux recovery errors (Section 2.2). 
This moderately boosts the error on the spectral index towards high 
S/N, but does not impact our source identifications. 

In addition to the bias from the steepness of the population, we 
also treat the fact that the source finder locates the maximum flux 
along RA and Dec., which provides an extra two degrees of free¬ 
dom. This can be corrected by finding the flux after two degrees of 
freedom are subtracted from the detection significance, following 
Vanderlinde et al. (2010). Note that this departs from the treatment 
in Crawford et al. (2010) and Vieira et al. (2010). Specifically, all 
sources are corrected by a factor (S/N)/ y/ S/N 2 — 2. For sources in 
the range 20-25 mjy, this is a correction of 0.5 percent at 148 GHz 
and 1.5 percent at 218 GHz. We have confirmed through simula¬ 
tions with synthetic sources implanted in the ACT data that this 
positional deboosting results in unbiased flux densities. Further¬ 
more, we have compared raw flux densities from the matched filter 
to flux densities from the ACT data at the positions of ATCA coun¬ 
terparts, when available. The latter should not be boosted due to 
maximizing the flux over position in RA and Dec. As expected, 
we find that positional deboosting accounts for the ratio between 


Table 2. Number counts, purity and completeness 2 . 


Flux range (Jy) 

N 

148 GHz 

Purity 

Completeness 

N 

218 GHz 

Purity 

Completeness 

Atync 

Adust 

0.015-0.02 

23 

100.0 ± 0.0 percent 

47.1 ± 3.2 percent 

8 

100.0 ± 0.0 percent 

10.7 ± 2.8 percent 

22 

0 

0.02-0.03 

47 

100.0 ± 0.0 percent 

75.3 ± 3.4 percent 

42 

92.9 ± 4.8 percent 

34.2 ± 3.5 percent 

46 

10 

0.03-0.05 

48 

100.0 ± 0.0 percent 

96.8 ±1.9 percent 

43 

97.7 ± 2.3 percent 

80.2 ± 3.0 percent 

48 

10 

0.05-0.09 

20 

100.0 ± 0.0 percent 

99.6 ± 0.6 percent 

21 

100.0 ± 0.0 percent 

98.7 ± 0.9 percent 

20 

1 

0.09-0.17 

13 

100.0 ± 0.0 percent 

100.0 ±0.0 percent 

12 

100.0 ± 0.0 percent 

100.0 ± 0.0 percent 

13 

2 

0.17-0.33 

4 

100.0 ± 0.0 percent 

100.0 ±0.0 percent 

2 

100.0 ± 0.0 percent 

100.0 ± 0.0 percent 

4 

0 

0.33-0.65 

2 

100.0 ± 0.0 percent 

100.0 ± 0.0 percent 

1 

100.0 ± 0.0 percent 

100.0 ± 0.0 percent 

2 

0 

0.65-1.39 

1 

100.0 ± 0.0 percent 

100.0 ± 0.0 percent 

1 

100.0 ± 0.0 percent 

100.0 ± 0.0 percent 

1 

0 

1.39-2.87 

1 

100.0 ± 0.0 percent 

100.0 ±0.0 percent 

1 

100.0 ± 0.0 percent 

100.0 ± 0.0 percent 

1 

0 


“The number of sources at each frequency with S/N = 5 over 455 deg 2 . Errors on number count are simply Poisson. See Fig. 6 for a graph of 
purity/completeness-corrected differential source counts. 

^Counts of synchrotron-dominated sources are taken relative to the 148 GHz flux density. 
c Counts of dust-dominated sources are binned according to their 218 GHz flux density. 
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the raw flux densities and the flux densities derived using ATCA 
counterpart locations. 

As a final consistency check, we note that for the sources observed 
by both ACT and Planck, flux densities are consistent at the ~1- 
2 per cent level, suggesting that the errors assumed here are very 
conservative and free from any systematic bias (Section 5.3). 


4.3 Purity and completeness 

The number of false detections at each frequency was estimated by 
running the detection algorithm on an inverted (negative tempera¬ 
ture) map in which we masked the sources and, in the case of the 
148 GHz map (for which the Sunyaev-Zel’dovich effect was non¬ 
zero), all ACT-detected and optically confirmed clusters of galaxies. 
With this approach, no spurious detections are found in the 148 GHz 
data down to an S/N of 5. In the 218 GHz data, four spurious de¬ 
tections at S/N < 6 were observed with raw flux densities in the 
25-31 mjy range (Table 2). 

The purity simulations appear to be consistent with our 
findings from source cross-identification (Section 5). At 
148 GHz, only one 5a source, ACT-S J023600-530237, is 
not cross-identified with other catalogues. At 218 GHz, five 
5 <t sources, ACT-S J024430-541605, ACT-S J035034-524801, 
ACT-S J062747—512614, ACT-S J063715-500414 and ACT-S 
J065207—551605, are not cross-identified with other catalogues, 
although we expect some of these to be real DSFGs, a hypothesis 
which will be tested through future follow-up observations. 

In order to estimate completeness of the catalogue for each band, 
we added 100 synthetic sources of a single flux density to the 
455 deg 2 region of the ACT maps used in this study. The map 
filtering and source detection algorithm (Section 3) were run, and 
the resulting catalogue checked for inclusion of the input sources. 
This procedure was repeated in each band in intervals of lOmJy in 
synthetic source flux density in the range 10-100 mjy. The com¬ 
pleteness at each flux density was interpolated to estimate the com¬ 
pleteness in flux density bins given in Table 2. 


5 COMPARISON TO OTHER CATALOGUES 

To further characterize ACT sources, we consider cross¬ 
identifications with other catalogues. Given the number of ACT- 
detected sources and the ACT beam size, around 1 out of 10 000 
sources randomly placed in the 455 deg 2 area considered here would 
coincide with an ACT source. Fig. 5 shows the ACT sources located 
in flux density space and cross-identifications from some overlap¬ 
ping catalogues with extragalactic source detections in the Southern 
hemisphere discussed in this section. The synchrotron-dominated 
sources occupy points near the power law S v oc v“ with spectral 
index a = —0.6. Dust-dominated sources will follow a line closer 
to a = 3.7, and are generally associated with sources detected with 
S/N >5 at 218 GHz only. 

As our cross-identifications come from radio catalogues at low 
frequencies, and dust-dominated sources lack strong synchrotron 
radiation, we expect fewer of the dust-dominated sources to be co¬ 
identified with a previously known radio source. Overlapping large- 
area, sensitive far-IR (FIR) catalogues would be the equivalent for 
our dust-dominated sources, but these do not exist. Even if they 
did, they would be too confused for positive cross-identification in 
this manner. At present, these unidentified sources are the subject 
of individual follow-up studies (e.g. Ikarashi et al. 2011; Scott et al. 
2011; Lupu et al. 2012; WeiB et al. 2013). 



Si 48 [ mJ y] 


Figure 5. The ACT sources cross-identified in flux density space. Sources 
have been selected with S/N > 5.0 in at least one frequency band (dotted 
lines). See the text for a more thorough description of cross-identification 
statistics with other catalogues. The diagonal lines follow the power law 
S v oc v a with a = 1.66 for the dashed line, distinguishing the a = —0.6 
(dash-dotted) synchrotron-dominated source population from the a = 3.7 
(dash-dot-dot) dust-dominated source population. 

5.1 General statistics of identifications 

Of our 191 sources, many have cross-identifications with sev¬ 
eral catalogues. 174 are identified with sources in the 0.84 GHz 
Sydney University Molonglo Sky Survey (SUMSS; Mauch et al. 
2003), and 122 of those cross-identified sources also belong to 
the 4.85 GHz Parkes-MIT-NRAO survey radio catalogue (Wright 
et al. 1994). 14 sources were cross-identified with the IRAS (Helou 
et al. 1988) at 12-100 pm. Sources were also identified with 
one or both of the AT20G survey (Murphy et al. 2010) and the 
1.4/2.0 mm SPT (Vieira et al. 2010) catalogues. Two of the re¬ 
maining eight unmatched sources were observed during the 2010 
November ATCA follow-up campaign (Section 2.3). This leaves 
six sources (ACT-S J023600-530237, ACT-S J024430-541605, 
ACT-S J035034—524801, ACT-S J062747-512614, ACT-S 
J063715—500414 and ACT-S J065207-551605) with no identi¬ 
fication in any of the available catalogues, all but one of which 
have dust-dominated spectra. The single unidentified source with 
a synchrotron spectrum, ACT-S J023600—530237, falls near our 
S/N = 5 threshold at 148 GHz and could be a false detection. 

A NASA/IPAC Extragalactic Database (NED) search using a 
radius of 2 arcmin cross-identified three sources from the ACT 
catalogue with galaxy clusters: Abell S0250, Abell 3391 and Sersic 
037/01. Increasing the search radius to 3 arcmin adds Abell 3128 
and Abell 3395 to this list. Simulations of the microwave sky suggest 
that only ~3 percent of galaxy clusters have their 148 GHz signal 
contaminated at 20 percent or more (Sehgal et al. 2010). 
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5.2 Matches to AT20G 

The AT20G survey, carried out from 2004 to 2008 at 20 GHz with 
a follow-up at 5 and 8 GHz, covered 6.1 sr of the southern sky to 
a flux density limit of 40mJy at 20 GHz (Murphy et al. 2010). Of 
our sources, 115 match sources in the AT20G catalogue, listed in 
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Table 4, all of which we classify as synchrotron dominated. Given 
the AT20G completeness limit and the predominance of flat spectra 
for our 148-GHz-selected sources (Section 7), faint ACT sources 
may not have matches in AT20G. We received time on the ATCA 
to measure flux densities for sources in the earlier Mil catalogue 
that did not appear in AT20G, the results of which are described in 
Sections 2.3 and 7. 

5.3 Matches to the Planck compact source catalogues 

The Planck satellite team has released two sets of all-sky catalogues 
in nine frequency channels, including two very close to those used 
by ACT. The first set of catalogues made up the Early Release 
Compact Source Catalog (Planck Collaboration XIII 2011). Very 
recently, a deeper catalogue based on 1.6 complete Planck sur¬ 
veys has been released, the Planck Catalog of Compact Sources 
(PCCS; Planck Collaboration XXVIII 2013). At 143 and 217 GHz, 
the Planck beam sizes are approximately 7 and 5, respectively. Like 
WMAP, Planck makes great gains by observing the entire sky simul¬ 
taneously at several frequencies, thus capturing all bright sources. 
On the other hand, due to their large beam sizes, these experiments 
detect only the brightest sources. In the case of the PCCS, the de¬ 
tection threshold at 143 and 217 GHz is ~ 400 mly; at lower flux 
densities the catalogue completeness drops well below 80 per cent. 
In the area treated in this study, Planck detects three of the brightest 
ACT sources. 

For a more complete comparison with the PCCS, we examined 
all bright ACT sources, combining the equatorial and southern re¬ 
gions mapped by ACT. Almost 50 matches were found in this 
wider comparison. At 148 GHz, the agreement between ACT flux 
densities and the PCCS flux densities, properly colour corrected 
and (slightly) extrapolated to match ACT’s central frequency, is 
excellent, at the 1-2 percent level, despite scatter introduced by 
source variability. At 218 GHz, source variability plays at least as 
large a role, and PCCS flux densities are ~5 per cent higher than 
ACT’s. If we remove a couple of the most variable sources, how¬ 
ever, the agreement of the 218 GHz flux density scales improves 
to & 1 per cent. At both frequencies, Planck fluxes are on average 
slightly higher. Agreement at the «l-2 percent level, however, is 
well within the expected uncertainty in the ACT flux densities (Sec¬ 
tion 2). This agreement suggests that our flux density values are free 
of systematic error. 

5.4 Matches to the SPT catalogue 

The SPT study of compact sources (Vieira et al. 2010) was based 
on observations of a square patch of sky of 87 deg 2 centred at 05 h 
right ascension, having only fractional ACT overlap. Nevertheless, 
2304 of the 3496 SPT candidate sources (those with S/N > 3 and 
flux densities >4.4 mjy) fall within the ACT survey region. We find 
32 cross-identifications with ACT sources when we search the SPT 

1.4 mm (220 GHz) and 2.0 mm (145 GHz) catalogues. 25 of these 
were 148 GHz selected and categorized in Vieira et al. (2010) as 
synchrotron dominated. The remaining seven were detected by ACT 
with S/N > 5 only at 218 GHz; three of the seven also have SUMSS 
or IRAS cross-identifications. To make a direct comparison of the 
flux densities of ACT-SPT cross-identified sources, we used our 
derived flux densities, but without introducing two angular degrees 
of freedom in the deboosting (as in Vanderlinde et al. 2010), as these 
were not incorporated by the SPT analysis (Crawford et al. 2010). 
We find an overall offset at the 7 ± 5 per cent level at 148 GHz and 
3 ± 7 per cent at 218 GHz, with ACT flux densities typically higher 


than those of SPT. At low flux, the ACT and SPT measurements 
of individual sources agree within their errors. However, at fluxes 
above 30 mly, individual sources may disagree at several sigma. We 
understand the discrepancy to be due to variability in radio sources, 
which is deferred to a future publication. 

Recently, using spectroscopic follow-up observations with the 
Atacama Large Millimeter/submillimeter Array (ALMA), Vieira 
et al. (2013) and Weifi et al. (2013) derived robust redshifts for 
about 18 of the SPT-detected dust-dominated sources from a larger 
sky area than in Vieira et al. (2010). One of these sources, found to be 
at redshift 5.66, matches ACT-S J034640—520505 which otherwise 
had no cross-identification. A second previously unmatched source, 

ACT-S J002707-500713, imaged with ALMA at 870 qm, is re¬ 
solved at the 0.5 arcsec scale, likely due to gravitational lensing, but 
awaits a spectroscopic follow-up. Four other ACT dust-dominated 
sources that were cross-identified with the Vieira et al. (2010) cat- a 

alogue now also have robust redshifts. ACT-S J055139—505800 | 

is at redshift 2.123, ACT-S J053250-504709 is at redshift 3.399, f 

ACT-S J052903—543650 is a source at redshift 3.369, with a lens o 

at redshift 0.140, and ACT-S J053817-503058 has a redshift of g> 

2.782, with a lens at a redshift of 0.404. It is likely that the few ACT B 

dust-dominated sources which lie outside the SPT footprint are also ,§ 

lensed DSFGs. 


Of the 24 ACT dust-dominated sources in the catalogue 
presented here, 18 of these are cross-identified with either 
IRAS (12 sources) and/or SPT (9 sources, including the re¬ 
cent ALMA observations). ACT-S J051506-534420, ACT-S 
J053311—523827 and ACT-S J055115-533435 were cross- 
identified with both catalogues. 10 of the dust-dominated sources 
are cross-identified with sources in the SUMSS catalogue, bring¬ 
ing the total for cross-identified sources to 19. The remaining 
five sources, ACT-S J024430-541605, ACT-S J035034-524801, 
ACT-S J062747—512614, ACT-S 1063715-500414 and ACT-S 
J065207—551605, have no matching counterpart, and signal-to- 
noise at 218 GHz in the range 5.27 < S/N < 6.35. Purity tests 
suggest that a couple of these could be spurious. However, given 
that three S/N ~ 6 sources are cross-identified with SPT, some of 
these sources are likely real detections, high-redshift galaxies lensed 


5.5 Matches to the IRAS catalogue 

Given the ACT beam size, a normal galaxy will be unresolved 
at redshifts z > 0.05 or distances greater than 200 Mpc. Conse¬ 
quently, only very nearby objects appear extended in our maps. The 
IRAS source population consists primarily of local (z <£. 1) dust- 
dominated ULIRGs; we thus expect our resolved sources to coin¬ 
cide with IRAS sources. We find 14 sources cross-identified with 
IRAS sources, out of a possible 829 IRAS sources that lie within 
the ACT survey area treated here. Two of these have synchrotron- 
dominated spectra, both of which are nearby galaxies. The 12 that 
show spectra dominated by dust re-emission are Galactic sources, in 
the Magellanic Clouds, or are known nearby star-forming galaxies. 
For example, for two local resolved galaxies particularly bright at 
ACT frequencies, NGC 1566 (ACT-S J041959-545622), a Seyfert 
two-arm spiral, and IC 1954 (ACT-S J033133-515352), a late- 
type spiral with a short central bar, ACT observes a higher 218 GHz 
flux density than at 148 or 20 GHz, confirming the dust-dominated 
nature of their spectra. 


5.6 Dust-dominated ACT sources 
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Figure 6. Differential number counts of ACT-selected sources. Derived from Table 2 and corrected for completeness, the ACT differential source counts are 
plotted together with models of radio and infrared source populations. The Planck Collaboration VII (2013) data points and the SPT data points of Vieira 
et al. (2010) are also plotted, both of which are consistent with ACT counts. The ACT and Planck points have Poissonian errors (la), whereas the SPT points 
include measurement and independent calibration errors. The FIR models on this scale predict number counts too low to be seen at 148 GHz. The ACT data 
are consistent with being dominated by radio sources at both frequencies. 


by intervening structure. Future follow-up observations are planned 
to clarify the nature of these detections. 

The main and supplementary 500 pm samples from the Herschel 
Multi-tiered Extragalactic Survey have lensing candidate densities 
of0.14 ± 0.04and0.31 ± 0.06 deg 2 , respectively (Wardlow et al. 
2013). An analogous sample from Herschel ’s H-ATLAS survey is 
0.35 ± 0.16 deg" 2 (Negrello et al. 2010). In the same bands as 
presented here, SPT finds 0.25 ± 0.02 deg -2 lensing candidates 
(Mocanu et al. 2013). This suggests that ~ 100 such sources exist 
in the field here, of which we see the high-flux tail. 

6 SOURCE COUNTS 

The completeness-corrected differential number counts for ACT 
sources based on the data in Table 2 are plotted in Fig. 6. For com¬ 
parison, completeness-corrected number counts from Vieira et al. 
(2010) and Planck Collaboration XIII (2011 ) are plotted as well. The 
ACT data fill in the flux density gap between the SPT and Planck 
catalogues at these frequencies, caused by the differences between 
experiments in sky coverage and sensitivity to point sources. The 
effect of calibration error is to shift the flux bins by ±2 percent 
(±2.4 percent) at 148 GHz (218 GHz; Section 2). 

The combined ACT, SPT and Planck total counts at 148 GHz 
and synchrotron-dominated source counts at 218 GHz are best fit¬ 
ted by the recently developed Tucci et al. (2011) C2Ex radio source 
model. This model divides the '“flat- or steep-spectrum blazar pop¬ 
ulation into BL Lac objects, with a region close to the AGN core 
dominating the observed emission (<0.3 pc), and flat-spectrum ra¬ 
dio quasars (FSRQs) with a break frequency indicating that the 
source of emission arises from an emitting region further from the 
AGN core (0.3-10 pc). The de Zotti et al. (2005) model for ra¬ 
dio source counts, while consistent with the counts below 0.1 Jy, 


overpredicts the counts at higher flux densities. Planck analysis 
finds that the de Zotti et al. (2005) model is consistent with their 
counts at frequencies up to 100 GHz, but overpredicts the counts at 
higher frequencies in the flux density region of «1 Jy, though they 
begin to suffer from incompleteness below Jy. Our few ^1 Jy 
brightest sources appear to be consistent with this finding, but lend 
little statistical significance. 

As well as the models for counts of radio sources, the source 
counts predicted for dusty starburst galaxies from Toffolatti et al. 
(1998), Lima et al. (2010), Bethermin et al. (2011) and Cai et al. 
(2013) are shown in Fig. 6. The brightest infrared sources in the 
Toffolatti et al. (1998) model are 10 mjy. Given that all sources in the 
ACT catalogue have flux densities greater than 10 mjy, these models 
predict that the 148- and 218-GHz-selected ACT catalogues should 
have few or no detected infrared sources. However, we find 24 
sources that have dust-dominated spectra with flux densities above 
20mJy, and Vieira et al. (2010) find 36 dust-dominated sources 
with flux densities above 15 mjy. While Toffolatti et al. (1998) do 
not incorporate the effects of lensing on observed number counts, 
the models of Lima et al. (2010), Bethermin et al. (2011) and Cai 
et al. (2013) modelled dust-dominated galaxies together with lensed 
high-redshift dust-dominated galaxies. However, the small number 
of dust-dominated sources present in Fig. 6 lends little constraining 
power, and we leave it to future larger area ACT studies to analyse 
the robustness of these models. The Bethermin et al. (2011) model, 
falling roughly in the middle of this group of models, was used for 
our flux deboosting. 

7 SOURCE SPECTRA 

Source SEDs can be used to differentiate source types by their dom¬ 
inant emission mechanisms. Assuming the commonly used simple 
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Table 3. Median spectral indices. 


Spectral index 

Synchrotron (all)" 

Synchrotron (Sws > 50 mjy) 

Synchrotron (Si48 < 50 mjy) 

Dust-dominated 

<*5-20 

—°.15ljj;^ 

-0.07ig^ 

-0.21+|g 

- 

<*20-148 

-0.42+°;^ 

—0.36«;“ 

-0.43+°;^ 

- 

<*148-218 

-0.55+°;“ 

-0.60ig-20 

— 0.51+qj® 

3 7+0.62 


"Quoted errors are the 68 percent confidence levels of the distribution. 


power-law model S{v) ocv“, a negative a is indicative of sources 
dominated by synchrotron emission, such as radio galaxies. Sources 
with free-free emission dominating will have an index close to 0. 
The high-redshift SMG population will have spectra dominated by 
re-emission of their prodigious optical and UV flux by the sur¬ 
rounding dust in a grey-body spectrum, with indices expected to be 
greater than 2 and more typically 3-4. 

We can divide the ACT source population according to sev¬ 
eral broad spectral groups: classical steep (and steepening) spec¬ 
trum sources, sources that peak within the frequency range under 
consideration, and sources that show flat, rising or upturned 
spectra. 

The ACT catalogue is dominated by synchrotron-dominated blaz¬ 
ers, which have variable flux densities. This variation in flux density 
is not periodic. For any single source, then, inferences about its spec¬ 
trum will depend on the epoch of observation, although not biased 
one way or the other. For the catalogue as an ensemble, however, a 
spectral study may give rise to insights about the average spectral 
behaviour of the galaxy populations. 

Table 3 summarizes the median spectral indices between pairs 
of frequencies for various subsets of the ACT data. In obtaining 
the spectral indices, we compare the deboosted flux densities from 
ACT with the raw flux densities from AT20G. The whole sample 
includes data from the 20 GHz ATCA follow-up observations taken 
in 2010 November (Section 2.3). 

Fig. 7 shows the 5-20, 20-148 and 148-218 GHz colour- 
magnitude diagrams. The radio-selected AGN (black points) are 
predominantly characterized by steepening of the spectra, with only 
a few characterized by extremely flat or inverted spectra. There is 
a clear trend towards more negative median spectral index with in¬ 
creasing frequency. Magenta points, showing rising spectra between 
148 and 218 GHz, denote sources ACT classified as spectrally dust 
dominated. They do not show up in the top two plots, indicating 
that these sources have flux densities falling below the detection 
threshold of the AT20G catalogue. 

The average spectral indices between AT20G and ACT frequen¬ 
cies indicate an underlying source population made up of FSRQs, 
a type of blazar, with AGN jet pointed along our line of sight (de 
Zotti et al. 2010). Ejected material flows through several shocked 
regions in the jet which locally enhance the radiation (Marscher 
& Gear 1985; Valtaoja et al. 1992). The observed spectral flatness 
is the superposition of many components with different turnover 
frequencies. At frequencies greater than approximately 100 GHz, 
however, Mil, Vieira et al. (2010) and Planck Collaboration XIII 
(2011) observe a steepening of the spectrum (ofi48_2i8 ~ —0.6) that 
until now had not been conclusively shown (Tucci et al. 2011). This 
change is possibly due to electron energy losses in the jet (‘elec¬ 
tron ageing’) or the transition to the optically thin regime in the 
extended radio lobes. The underlying physical mechanisms have 
been contested for more than a decade and remain the subject of 
ongoing study (e.g. Ghisellini & Tavecchio 2008; Nieppola et al. 
2008; Sambruna et al. 2010). 



Figure 7. Colour-magnitude diagrams comparing the 5-20 GHz (top), 20- 
148 GHz (middle) and 148-218 GHz (bottom) spectral indices for ACT- 
ATCA cross-identified sources. The synchrotron-dominated radio galaxy 
population is dominated by sources which have consistently falling SEDs 
towards higher frequencies. 

All but one of the synchrotron-dominated ACT detections with 
flux density >50mJy have cross-identifications in AT20G. Below 
this flux density, the mean spectral indices of the population of ACT- 
AT20G cross-identified sources are biased towards the negative by 
the incompleteness of the AT20G catalogue below 100 mjy at 5 GHz 
and below 40 mjy at 20 GHz. It is therefore illustrative to further 
divide this subpopulation according to the 148 GHz flux density. 

Fig. 8, a radio colour-magnitude diagram, plots the spectral in¬ 
dex 0(20—148 against the 148 GHz source flux density. Black and 
grey points identify the Ml 1 sources that were cross-identified with 
the AT20G catalogue. Black points are for sources with flux densi¬ 
ties >50 mjy, which represented a complete sample. Grey points de¬ 
note the fainter, incomplete sample. The blue points were obtained 
by calculating the spectral index for the previously unmatched, 
lower flux density subsample, followed up with ATCA (see Table 4 
for AT2G0 source IDs with an asterisk). The population represented 
by the blue points fills in the picture remarkably for the fainter (be¬ 
low the 50 mjy flux density at 148 GHz, below 40 mjy at 20 GHz) 
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Figure 8. Radio colour-magnitude diagram using 20-148 GHz spectral in¬ 
dices for ACT-AT20G cross-identified sources. The data are divided between 
flux densities at 50 mjy at 148 GHz. The low-flux sample w; n implete 
and suffered from selection bias that favoured sources with more nega¬ 
tive spectral indices. Data from the ATCA follow-up study for the Mil 
ACT 148 GHz sources not cross-identified with AT20G are shown as green 
squares. The high-flux sample, denoted with black points, has a median 
spectral index —0.37. Prior to ATCA follow-up observations, the lower flux 
ACT sources cross-identified with the AT20G catalogue (grey crosses) had 
a median spectral index of —0.52. Including the follow-up data, the more 
complete lower flux sample has a median spectral index of —0.43. 


population, especially in the region with spectral index of approxi¬ 
mately zero. 

The unbiased S1 4 g > 50mJy sample has a 20-148 GHz median 
spectral index of —0.36iJ5;|J. For the 514s < 50mJy subsample 
prior to a follow-up with ATCA (grey points only), the 20-148 GHz 
spectral index was <*20-148 — — 0 . 52 +Q^g. However, with the addi¬ 
tion of the synchrotron-classified ACT-selected sources with 2010 
ATCA follow-up data (grey and blue points), the full sample has 
20-148 GHz spectral index 0*20-148 = —0.42+^. This supports the 
hypothesis that the lower flux contingent is probing the same pop¬ 
ulation of synchrotron-dominated sources (blazers), and that there 
is not much evolution of their spectral index with flux. The full 
data set suggests that the increased scatter at lower fluxes is due to 
variability (which will have a larger relative effect on the fluxes) 
and decreased S/N from flux error. 

For our dust-dominated sources, we find a median spectral index 
of ai48-2i8 — 3-7^a86- F° r sources detected above 5er at both 150 
and 220 GHz, SPT derives £*150-220 = 3.35 ± 0.7 (Mocanu et al. 
2013). These indices are consistent with the findings of Dunkley 
et al. (201 1), where the effective index of unresolved IR emission 
as determined from ACT data is 3.69 ± 0.14, and the best-fitting 
mean spectral indices of <*f 5 0 _ 22 o — 3.86 ± 0.23 for Poisson- 
distributed sources and <*f 50 _22o = 3-80 ± 1.3 for the clustered 
component, found below the SPT detection threshold (Hall et al. 
2010), a population expected to be dominated by dust-dominated 
sources. If we model dust emission as a modified blackbody, in the 
Rayleigh-Jeans (RJ) limit, the dust emissivity index is related to the 
effective spectral index as fi = a — 1 — 1.7, consistent with models 
of ~30 K dust made from graphite and silicate grains (Draine & 
Lee 1984). We leave a more rigorous analysis involving redshifted 


grey bodies, where the RJ limit is not as good an approximation, 
and the implications for star formation to future work. 


8 CONCLUSION 

We have described the extragalactic source population at 148 and 
218 GHz found in a 455 deg 2 region of the ACT 2008 southern 
strip , centred on declination —52?5. This updates the ACT 148 GHz 
catalogue by using the data released with D13, extends the results 
in Ml 1 to two bands and treats noise as more local, which in turn 
yields a higher S/N. We detect 191 sources above S/N = 5 in at 
least one of the ACT 148 and 218 GHz frequency bands, spanning 
flux densities 14—1700 mjy (Table 4). Known redshifts are as high 
as ~6, with measurements ongoing. The catalogue is estimated to 
be 100 percent pure and 96.8 percent complete above 30mJy at 
148 GHz, and 97.7 percent pure and 80.2 percent complete above 
30 mjy at 218 GHz. We have confirmed flux recovery of the pipeline 
and jointly deboosted the flux densities in both bands. 

The multifrequency nature of our observations allows for in¬ 
ternal classification of sources into three broad classes of sources 
based on their spectra: synchrotron-dominated sources (the vast 
majority of which are cross-identified with radio catalogues), low- 
redshift dust-dominated sources with IRAS counterparts (typically 
ULIRGs) and dust-dominated sources also observed by SPT or with 
no cross-identification, either shown or assumed to be high-redshift 
star-forming galaxies. This last class, only recently observed at mil¬ 
limetre wavelengths, has many of the properties expected of the 
progenitor population of massive, modern-day, elliptical galaxies, 
background SMGs whose flux has been magnified through gravita¬ 
tional lensing by a foreground galaxy or galaxy group. This inter¬ 
pretation, bolstered by population synthesis analyses (e.g. Thomas 
et al. 2005), is being validated with follow-up observations. 

A comparison with other catalogues shows that 97 percent of 
ACT-detected sources correspond to sources detected at lower or 
higher frequencies. The 148 GHz source counts are fitted reason¬ 
ably well by the C2Ex radio model of Tucci et al. (201 1), the most 
current model for radio sources. According to the analysis of the av¬ 
erage spectral indices derived from the combined AT20G and ACT 
data sets, the ACT data support the case for a spectral steepening 
towards higher frequencies above 100 GHz for AGN. The ACT dust- 
dominated source population has a median spectral index, <*148-218, 
of 3.7 ^q gg. Properly linking these sources into the broader context 
of galaxy formation and evolution is of cosmological interest, and 
a goal of future work. 

The analysis presented here uses only the 2008 data of ACT’s 
southern strip, representing only l/6th of the data ACT obtained 
between 2007 and 2010. In future work, we will extend our analysis 
of the source population to include the full data set integrating both 
the southern and equatorial regions observed by ACT. The ACT 
equatorial strip overlaps with deep Sloan Digital Sky Survey Stripe 
82 observations (Annis et al. 2011). Thus, as well as increasing the 
sky coverage and number counts for the ACT sources, future work 
will include joint analyses with optical data. 
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